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Abstract 

A multiple-scale eddy viscosity model is described in this 
paper. This model splits the energy spectrum into a high 
wave number regime and a low wave number regime. Divid- 
ing the energy spectrum into multiple regimes simplistically 
emulates the cascade of energy through the turbulence spec- 
trum. The constraints on the model coefficients are deter- 
mined by examining decaying turbulence and homogeneous 
turbulence. A direct link between the partitioned energies 
and the energy transfer process is established through the 
coefficients. This new model has been calibrated and tested 
for boundary-free turbulent shear flows. Calculations of 
mean and turbulent properties show good agreement with 
experimental data for two mixing layers, a plane jet and a 
round jet. 


1 Introduction 


In turbulent flows, the mean flow performs deformation 
work which tansfers energy from the mean flow to the large- 
scale turbulent eddies. The turbulent kinetic energy con- 
tained within these large eddies is passed to the smaller ed- 
dies by vortex stretching. Once the energy has been passed 
into eddies near the Kolmogorov scale, it is then dissipated 
by the molecular viscosity. This process can be thought of 
as a turbulent kinetic energy cascade [1], In other words, the 
turbulent kinetic energy is passed through the wave n um ber 
spectrum as it cascades from large to small eddies. 

The multiple-scale turbulence model which is proposed in 
this study splits the energy spectrum into low and high wave 
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number regions. The low wave number region contains the 
large-scale eddies and the high wave number region con- 
sists of the smaller, less energetic eddies. A division of this 
nature models the cascade of energy from the production 
region where the energy is initially created by the mean 
straining work to the dissipation region where thin energy 
is eventually dissipated. This concept is illustrated in Fig- 
ure 1. Although a simple two-part division of the energy 
spectrum cannot fully model the cascade of energy, it can 
simulate the nonequilibrium energy transfer process which 
is beyond the capability of all single-scale models. 

This concept was incorporated into an earlier multiple-time 
scale turbulence model by Hanj&lic, Launder and Schiestel 
[2] (hereafter HL&S). Based on the same modeling method- 
ology that was used in the development of the “standard” 
k-t model [3], HL&S derived four transport equations to 
desribe the turbulent characteristics of the two regions. 
Consequently, HLicS formulated two transport equations 
for the partitioned turbulent kinetic energies. They also de- 
veloped transport equations for the rate of energy transfer 
between the two scales and for the rate of energy transfer to 
the small-scale eddies. In this model, the coefficients were 
written as functions of both the ratio of the partitioned en- 
ergies and the ratio of the spectral energy transfer rates. 
A term containing the mean vorticity was included in the 
energy transfer rate equation to account for the increased 
energy transfer rates in irrotational shear flows. Their re- 
sults [2] showed fairly good agreement between the model 
predictions and the experimental data for jets and boundary 
layers. 

Kim and Chen [4] (hereafter K&C) developed another 
multiple-scale model based on the energy partitioning idea 
introduced by HL&S. In this model, the transport equa- 
tions were modified to include an extra source term in both 
the energy transfer rate equation and the energy dissipation 
rate equation. In addition, the turbulent velocity scale was 
characterized by the total turbulent kinetic energy. K&C 
calibrated the model constants which appear in the energy 
transfer equations for simple, wall-bounded turbulent flow 
problems. This model has been used for several boundary- 
layer flow problems using the appropriate near-wall correc- 
tions. 
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The model in this study is also based on the energy par- 
titioning concept of HL&S. In particular, the model coeffi- 
cients are dynamically dependent upon the partitioning of 
the energy spectrum. The variable nature of the coefficients 
adjusts the model to different flow situations. There is no 
need for the extra source terms used in K&C’s model or the 
rotational straining term in the HL&S model. These model 
coefficients have been calibrated for homogeneous shear flow 
and decaying grid turbulence. The present multiple-scale 
turbulence model has been tested for boundary-free shear 
flows. Two mixing layers at different speed ratios, a planar 
jet and a round jet have been evaluated. All the computa- 
tions show reasonably good agreement with the data. 


2 Model Equations 


Mean Flow Equations 


For incompressible turbulent flow, the ensembled-averaged 
equations for continuity of mass and momentum are written 
as 



(i) 


DUi d ( dUi \ 1 sp 

Dt dzj \ dzj * ) pdzi 


( 2 ) 


where — ' ujuj 18 the turbulent Reynolds stress tensor. Using 
the eddy viscosity concept, the Reynolds stress can be re- 
lated to the mean strain rate and a turbulent eddy viscosity, 


- U{Uj = V t 


(dn.duA 

\dZj + dZi ) 



Now, the momentum equation can be written as 



(3) 

(4) 


The turbulent eddy viscosity, v ti can be characterized by 
the local turbulent kinetic energy and the local length scale 
of the energy containing eddies, 

v t cckll. (5) 

The definition used for this length scale is the primary dis- 
criminating factor of eddy viscosity turbulence models. For 
instance, in a k-t model this length scale is written in terms 
of the turbulent kinetic energy, Jt, and its dissipation rate, e, 


In the present model, the length scale is described by the 
total turbulent kinetic energy and the rate of spectral en- 
ergy transfer, i.e. / = The turbulent kinetic energies 
and the rates of energy transfer are determined by modeled 
transport equations similar to the "standard” k-t equations. 
These transport equations are described in the following sec- 
tion. 

Turbulence Equations 

In this multiple-scale turbulence model, the energy spec- 
trum has been split into a region where the turbulent kinetic 
energy has been produced by interaction with the mean flow 
and a region where the turbulent energy has been trans- 
ferred from the production region. This division can be 
graphically represented by Figure 1. Now, kp is the kinetic 
energy contained within the production region and c p is the 
rate at which energy is passed from the low wave number 
range into the high wave number range. At the high wave 
number end of the spectrum, kt is the kinetic energy con- 
tained in the smaller eddies and 6| is taken to be equivalent 
to the dissipation rate. The modeled transport equations 
are 



Here, Pu p is the production of kinetic energy by the large 
scale eddies. Note that the e p term serves as a sink in the 
equation for kp and a source in the kt equation. The source 
and «jnlc terms in the energy transfer rate equations (c p 
and tt) are related to those in the turbulent kinetic energy 
equations by their corresponding time scales, i.e. ^ and 

These transport equations, equations (7) through (10), 
posses several differences from the earlier multiple-scale 
models. HL&S include a rotational straining term in the 
€ p energy transfer equation (8) to improve their model’s 
performance in axisymmetric flows. Based on dimensional 
r easoning , K&C [4] include additional terms in the c p and 
ct equations which are nonlinear in Pkp and e p , respectively. 
The present model is the least complicated as it uses neither 
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the rotational straining term nor the nonlinear production 
terms for energy transfer rates in the model equations* 


Accounting for the multiple partitions in the turbulent ki- 
netic energy spectrum and the spectral energy transfer rates, 
the eddy viscosity in the present model is defined as 


t/t = c M 


(ti + k p ) 2 

€ r 


( 12 ) 


and Cp = 0.09. This is the same relationship for eddy vis- 
cosity used by KfcC. With this formulation, the single scale 
eddy viscosity model will be recovered when c p approaches 
€ in equilibrium turbulence. The <r coefficients in equations 
(7) through (10) are assumed to be the following constants, 


= aki = 1*0 and = <r tt = 1.3. (13) 

The other coefficients in equations (7) through (10), namely 
Cpi, Cp 2 , c n and c^, are the modeling coefficients discussed 
in the following section. 


3 Model Coefficients 


The coefficients for this incompressible model have been de- 
termined from analyses of homogeneous and decaying tur- 
bulence. 

Grid Turbulence 


In homogeneous decaying grid turbulence, the turbulent 
quantities are functions of time only and equations (7) 
through (10) can be simplified to 


dk p 

dt 

& 

1 

II 

(14) 

de p 

dt 

= -*4 

Kp 

(15) 

dk t 

dt 

vCT 

1 

& 

II 

(16) 

dc t 

dt 

- r c « £ p - «* 

— cti— ctar-. 

k t k t 

(17) 


Most of the experimental evidence suggests that the turbu- 
lent kinetic energy decays in time and can be represented 


by 

kpo \W 

(18) 

and 

ijir 

kto \to J 

(19) 


where n is the decay rate and is typically of the order 1.2. 
From the above kinetic energy equations, (14) and (16), the 
energy transfer rates must decay as 


ilL 

€p0 



( 20 ) 


and 



( 21 ) 


Therefore, the Cp* coefficient from equation (15) can be re- 
lated to the decay rate, 


Cp 2 = 


n + 1 
n 


( 22 ) 


Manipulating the simplified transport equations, (14) 
through (17), and using the relations given by equations 
(18) through (21) yields a relationship between c % 2 and the 
other coefficients, 


C 12 = 


*1 + &P3 

fc + 1 


(23) 


Homogeneous Shear Flow 

Guidelines e*n be established for determining the remain- 
ing two constants, Cpi and C\ 1 , by examining the physical 
behavior of homogeneous shear flows. In this flow situation, 


the turbulent transport equations reduce to 



dkj, _ _ 

dt ~ Pkp ~ p 

(24) 


de p c p £p 

~ e * i k; p *~ c »k. 

(25) 


dkt 

It = ep - et 

(26) 


dt, t,i, c? 

— = C—-C-. 

(27) 

If we define 

ft 

II 

(28) 

and 

M '9 
11 

(29) 

then, equation (24) can be combined with (25) to give 



c v dkp (a - 0) 

k p dc p CpiQ — Cpi0 

(30) 


Likewise, equations (26) and (27) can be combined with the 
following result, 

= £zi— . (31) 

k t de t ctifi — Cxi 
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Assuming that P and the percentages of the kinetic energy 
contained in k p and k t remain nearly constant, then 

£p dk, _ ft * _ ( a ~ft)_ (32) 

j fe p dc p * de , cpia - <* 2 /? 

= (33) 

i, de, i de, cnp - c« 


and 


the model has been compared to experimental data 

and to the multiple-scale models ^Hanj^c, L«m « « 

Schiestel f2l (HL&S) end K™ 80(1 01101 M ( K&0 ) . 
tte Candid” k-e model. In all cases the parabohe ^sota- 
tne bww» ^th an initial plane and the flow 

S^» th^P^io“ margin the axial dire. 

i^lutions are decked to insure that they mamtain 
asdf-perserving profile. The results of : thesec Rations 
presented in Figures 2 through 18 and m Table 2. 


From the experiment of [5], the ratio ££ has been 

l Ctorly.thi. 4.1 Planar Mixing Layers 


to be 1 085. In a similar experiment, xavouiu™ 

[61 found this ratio to be between 

ratio is on the order of one. For simplicity, this term, T 3J 7 
, has been assumed to be umty. Th«efore, the foUoi^g 
expressions for c,i and ci can be found from equations (32) 

and (33), 


Cpl 


«HK- 


_ • • _ lavnra tHft flow is tO hsVC 

air and the jet. ift o gwu flow 

turbulence quantities are applied at the edges oi 
field and an equally spaced grid is used. 


_ ft -1 . *1 

^ ~ p + P ' 


(34) 

(35) Speed ratio = 0.0 


Equation (23) is now, 


0 - 1 + cpaftit 
Ct5= ft + ftfe - 1 


(36) 


and, Cp 2 is defined in equation (22). The coefficients, Cpr.Ci 

80 . .• „ ft* it and r*-. In the present 

and e, a are functions of u , «, “ 17 

model, the ratios described in equations (28) and (29) are 

assumed to be the following constants, 

(37) 


a= ?h 1 = 2.2 

e< 


3 = 5t = 1 . 05 . (38) 

e< 

These constants have been calibrated considering that ex- 
perimental measurments of homogeneous shear flow suggest 
that the ratio f should be near two. These coefficients are 
summarized in Table 1. 


Notice that the tone of i. to v»y - to «*» 

of the turbulent kinetic energy m the small scales to the 
energy contained in the large scales changes Since most of 
the energy is contained in the large scales, this ratio ahould 
£. than one [1]. This coefficient adaptabih^r al- 
lows the model to adjust to different flow situations and ^ a 
unique characteristic of multiple-scale eddy viscosity mod- 

els. 

4 Results and Discussions 

The present multipl.scale model has been tested for tiro 
planar mixing layers and two jet flows. As part of these 


Tn Fisure 2 the three multiple-scale models and the “stan- 
? model are compared to the experimental data 

i y ^o“ d of" “ dMu ^ve 

than eitherPatel’s data or the computations. The sp 
than eltb , , th four ^dy viscosity models and 

ing rates predicted by thetour y 
the spreading rate measured by Patel areJ hstea ^ 

For mi-ring layers, the spreading rate is defined " <i* “ 

“standard” k-t model are closest to the data. K&t, s moa 
significantly underestimates the growth rate of the mixing 

layer. 

Th. tendency of K&C-. model to -£*^*££5 
f: cure- and to tom 

tom ett«-. Apparently, to cumjmton^ 

ZZ. toft further tom.de to lo. ep««d .de of to 

£T«Sto experimentally toto — >*»»• 


t‘ 7^ ttoS to mitog layer -her. to turbulence 

^ectrum i. in equilibnum. The tot that {j remrin. mui 
L thm, one indtote. that mo* of the enmgy » eo.tm.rf 
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by the large scale eddies as stated by Tennekes and Lumley 

[!]• 

Speed ratio = 0.3 

The next case considered is a mixing layer with a speed ratio 
of 0.3. In Figure 6, the mean velocity profiles predicted by 
all four models are in very good with the experimental data 
measured by Spencer et. al. [9]. The present model and 
the "standard” k-e model predict an almost identical ve- 
locity profile. HL&S’s and K&C’s models under-predict the 
growth of the mixing region which is indicated by the nar- 
rower mean velocity profile and also by the lower spreading 
rates listed in Table 2. 

Looking at turbulent quantities, the "standard” k—e and the 
present k p -t p -k t -tt model predict the kinetic energy distri- 
bution very well as can be seen in Figure 7. Although the 
models due to HL&S and KfcC yield good mean velocity 
predictions, they both under-predict the peak turbulent ki- 
netic energy level. The present model and the "standard” 
it — e model yield very good predictions of the turbulent ki- 
netic energy profile. All the eddy viscosity models in this 
study under-predict the peak shear stress shown in Figure 
8. Here, the "standard” k-t and the present model do the 
best job of predicting the shear stress profile however, they 
are both slightly low in their predictions of the peak value. 
Again, the other two models are quite low in their predic- 
tions of the peak turbulent shear stress. 

Grid Resolution Analysis 

Since the mixing layer with a speed ratio of 0.3 gives the 
thinnest shear layer and the slowest growth rate, this case 
would be the most dependent upon grid resolution. For 
all the calculations discussed thus far there were 65 points 
across the flow field. Figures 9 and 10 compare the solutions 
obtained with 65 points to the same calculation with 101 
points across the domain. As these Figures indicate, the 
results are essentially unchanged as the grid is refined. 


4.2 Planar Jet 

For the jet flow simulations, the initial plane is split with a 
uniform velocity and kinetic energy profile comprising ap- 
proximately half of the domain and quiescent air compris- 
ing the other half. Zero gradient boundary conditions are 
applied at the centerline of the jet and the grid has been 
clustered towards the centerline to improve the accuracy of 
this boundary condition. 

As can be seen in Figure 11, the mean velocity profiles pre- 
dicted by the present model, K&C’s model and the "stan- 


dard” k-t model lie within the scatter of the experimental 
data of [10], [11] and [12]. The present multiple-scale model 
is particularly close to the data of Heskestad [12]. Again, 
the HL&S model under-predicts the turbulent growth rate 
which can be inferred by the narrow mean velocity profile 
and the predicted spreading rate. (During this study, the 
HL&S model has been found difficult to use. The model ap- 
pears to be sensitive to initial conditions, marching step rise, 
etc.. This is especially true for the planar and axisymmet- 
ric jets flows.) With the exception of HL&S’s model, there 
is very little disagreement between the predicted spreading 
rates and experimental spreading rate given in Table 2. For 
jet flows, the spreading rate is defined as the rate of change 
of the half velocity point, i.e. -J-j 1 . The predicted spread- 
ing rates range from 0.103 by K&C’s model to 0.114 by 
the "standard” k-t model. All these spreading rates are in 
reasonably good agreement with the experimentally deter- 
mined rates. 

There is considerably more disagreement between the mod- 
els for the predicted turbulent kinetic energy levels, shown 
in Figure 12. HL&S’s model under-predicts the kinetic en- 
ergy profile which is not surprising given the low predicted 
spreading rate. Both the present model and K&C’s model 
predict relatively flat profiles which match the data of Hes- 
kestad [12]. The "standard” k-t model is in very good agree- 
ment with the experimental data, being especially close to 
Bradbury’s measurements. 

Both multiple-scale models and the "standard” k—t model 
correctly predict the peak shear stress shown in Figure 13 
although they tend to exaggerate the width of the jet. The 
model by HL&S under-predicts the peak shear stress, again 
due to the under-prediction of the growth rate for the jet. 

Figure 14 shows the ratios of and ^ through the jet for 
the present multiple-scale model. Notice that these ratios 
are nearly constant and are therefore consistent with the 
assumption made in deriving the modeling constants. The 
ratio, f*-, is near one, except near the centerline of the jet, 
indicating that the energy spectrum is in equilibrium. On 
the other hand, the ratios of kinetic energy are much smaller, 
indicating that most of the energy is in the larger scales. 


4.3 Round Jet 

Using the "standard” k-t turbulence model often gives poor 
predictions for the spreading rate and flow properties of a 
round jet. There are several methods to adjust the “stan- 
dard” k—t model, e.g. by chang in g the coefficient in the 
dissipation equation or by adding a vortex stretching term 
[13]. The "standard” k-t model predictions given in Fig- 
ures 15 through 17 have applied no correction methods. As 
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previously mentioned, a rotational straining term is used scales crudely emulates the cascade of energy from its pro- 
by WTJb-S in their multiple-scale model to decrease the en- duction by the mean flow to its eventual dissipation by the 
ergy transfer rate from the production region in rotational molecular viscosity. This modeling concept can simulate the 
flows. Neither the present model nor K&C’s model uses a nonequilibrium energy transfer process. 


three-dimensional correction term; yet, as shown in Figure 
15, they give the best match to the experimental data of 
Rodi [14] and Wygnanski & Fielder [15] for mean velocity. 

The turbulent kinetic energy is over-predicted by the “stan- 
dard” k-€ model and under-predicted by HL&S’s model 
as shown in Figure 16. Closer to the data, but still not 
giving good predictions, are the present model and K&C’s 
model. Both models predict the correct centerline turbu- 
lent kinetic energy, however, the outer edges of the jet are 
over-predicted. None of the models in this report does a 
good job of predicting the turbulent kinetic energy for an 
axisymmetric jet. 

Figure 17 shows that there is also a wide spread in the shear 
stress predictions. Without adjusting the modeling coeffi- 
cients or adding correction terms the “standard” k-t model 
predicts shear stress levels almost double those given by 
the data. Both the present model and K&C’s model are 
relatively dose to the data compared to the “standard” k- 
€ model. Although, both models over-predict the peak level 
by approximately 25 percent. HL&S’s model predicts the 
narrowest shear stress profile in Figure 17 compared to the 
other models. Not surprisingly then, the HL&S model gives 
the best prediction for the spreading rate. The other two 
multiple-scale models, i.e. the present model and K&C’s 
model, predict spreading rates which are 17 and 24 percent 
too large, respectively, in comparison to the data. 

Looking at the ratios of ^ and shown in Figure 9, the 
flow field is not in equilibrium and the energy transfer to 
the dissipative scales increases near the centerline of the jet. 
The “standard” k-e model has no means to account for the 
increase in energy transfer rate near the centerline of the jet 
and consequently over-predicts the turbulent kinetic energy 
and the spreading rate. 


5 Conclusions 

A multiple-scale eddy viscosity turbulence model with solu- 
tion dependent coefficients has been developed. This new 
model splits the energy spectrum in two regions thereby in- 
troducing two characteristic length scales. One length scale 
is associated with the turbulent large-scale eddies which are 
responsible for generating turbulent kinetic energy by in- 
teractions with the mean flow while the other length scale 
corresponds to the smaller and less energetic eddies. This 
partitioning of energy between the large scales and the small 


In this paper, the multiple-scale turbulence model has been 
calibrated wring decaying grid turbulence and homogeneous 
shear turbulence. This model has been tested against two 
miTing layers at different speed ratios, a planar jet and an 
axisymmetric jet. For all cases, the mean flow and turbulent 
properties have compared reas on ably well with the experi- 
mental data. 
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Cp - 0.09 

n = 1.2 

a =2.2 

P = 1.05 

Cpl = (l - f ) + f 

C P J - n 

*1 = ^ + 7 

0— l+e*30rj- 
0,3 ~ 1 

*k r = 1.0 

<Tk, = 1.0 

&t r = 1-3 

<r tf = 1.3 


Table 1: Multiple-Scale Turbulence Model Coefficients 



Experiment 

HL & S 

kp-€ p -ki~€ t 

K&C 

kp—Cp-kt-tt 

present 

standard 

k-t 

Planar Jet 

0.11-0.12 

0.076 

0.103 

0.104 

0.114 

Round Jet 

0.085-0.095 

0.106 

0.118 

0.111 

0.126 

Planar Mixing Layer, r=0.0 

0.179 

0.146 

0.126 

0.152 

0.159 

Planar Mixing Layer, r=0.3 

0.052 

0.061 

0.064 

0.078 

0.082 


Table 2: Spreading Rate Comparisons for Free Shear Flow 
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Figure 1: Energy spectrum splitting for the multiple-scale Figure 3: Turbulent kinetic energy profile for a planar mix- 
mcK j e j ing layer with a speed ratio of 0. Umax'- maximum velocity. 



Figure 2: Mean velocity profile for a planar mixing layer Figure 4: Shear stress profile for a planar mixing layer with 
with a speed ratio of 0. Umax • maximum velocity. a speed ratio of 0. Umax- maximum velocity. 
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Figure 5: Energy transfer rate and kinetic energy ratios for Figure 7: Turbulent kinetic energy profile for a planar mix 
a plane mining layer with a speed ratio of 0. mg layer with a speed ratio of .3. AU = Umax — U m in • 
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Figure 6: Mean velocity profile for a planar mixing layer Figure 8: Shear stress profile for a planar mixing layer with 
with a speed ratio of .3. Umax ■ maximum velocity. a speed ratio of .3 . AU = U max — U m in- 
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Figure 9: Grid resolution analysis of the mean velocity pro- Figure 11; Mean yeiodty profile for a turbulent, planar jet. 
file for a planar mixing layer with a speed ratio of .3. Umax : . centerline velocity, 

maximum velocity. 



Figure 10: Grid resolution analysis of the shear stress profile Figure n . -jn^ent kinetic energy profile for a planar jet. 
for a planar mixing layer with a speed ratio of .3. AU = „ centerline velocity. 

Umax ~ Umin - 
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Figure 13: Shear stress profile for a planar jet. Umax • cen- 
terline velocity. 


Figure 15: Mean velocity profile for a round jet. Umax ’ 
centerline velocity. 



Figure 14: Dissipation and kinetic energy ratios for a planar 
jet 



Figure 16: Turbulent kinetic energy profile for a round jet 
Umax - centerline velocity. 


11 



l^/kp ond * ( /c p i»v«100/Umox**2 



y/* 

Figure 17: Shear stress profile for a round jet. Umat v* cen- 
terline velocity. 



Figure 18: Dissipation and kinetic energy ratios for a round 
jet. 
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